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Abstract: In this paper we formulate a corporate bond (CB) pricing model 
for deriving the term structure of default probabilities (TSDP) and the recov- 
ery rate (RR) for each pair of industry factor and credit rating grade, and 
these derived TSDP and RR are regarded as what investors imply in forming 
CB prices in the market at each time. A unique feature of this formulation is 
that the model allows each firm to run several business lines corresponding to 
some industry categories, which is typical in reality. In fact, treating all the 
cross-sectional CB prices simultaneously under a credit correlation structure 
at each time makes it possible to sort out the overlapping business lines of the 
firms which issued CBs and to extract the TSDPs for each pair of individual 
industry factor and rating grade together with the RRs. The result is applied 
to a valuation of CDS (credit default swap) and a loan portfolio management 
in banking business. 



1. Introduction 

In the financial industry credit risk is a central theme of interest in banking, in- 
vestment and derivatives. In the past financial economic theory and no-arbitrage 
theory in mathematical finance has provided many important results and thoughts 
for looking into many credit-risk-related problems. These normative theories are 
very important as a benchmark for those in financial industry. However, a blind 
application of the theories to the real world may cause a big problem, which hap- 
pened in the 2008 Financial Crisis. In fact, rating agencies wrongly applied the 
theories to valuation of mortgage-backed securities and expanded the serious Sub- 
prime Problem, while insurance companies issued credit risk swaps beyond their 
capacity of possibly taking risks and made a serious impact on the US economy, 
which in turn has been badly affecting world economies. In this paper, the author 
often emphasizes the importance of testing the empirical validity of "theories" and 
the limitation of the framework of no-arbitrage theory from a practical viewpoint 
below. But he never intends to deny the importance and value of theory itself, 
and thinks it important to now discuss openly about what caused the crisis among 
academics and practitioners in order not to let it happen again. 
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In credit risk analysis in mathematical finance, a time-continuous model is usu- 
ally formulated together with no-arbitrage concept and it is often the case that 
a spot rate process is used for instantaneous interest rates, the default-event gen- 
eration process is associated with such a model as hazard model or Merton-type 
model and the instantaneous recovery rate at each moment is constant. The result- 
ing model is usually Markovian and univariately stochastic, it individually evaluates 
the default probability of each product or each firm and the data used in applica- 
tions is mostly past default rates, rating scores, and firms financial data together 
with stock prices. It should be pointed out that the Markovian properties, which 
are the results of the assumption of time-continuous diffusion model, do not hold 
for credit risk processes in reality in general. It is because business cycles that 
are closely associated with credit variations of firms are not Markovian, which is 
a limitation of the model. On the other hand, Kariya and Liu (2004) treated a 
time-discrete theory to include non-Markovian models. 

In addition, since credit risk in reality is related to many factors such as indus- 
try business cycles which may depend on import-export relations, exchange rates, 
resource prices, etc, it does not admit a complete model except for the case in 
its theoretical assumption, meaning that the no-arbitrage concept in mathematical 
finance is not necessarily effective in reality. This in turn implies that in practice 
there is no unique valuation via the no-arbitrage. 

In this paper we formulate a corporate bond (CB hereafter) pricing model that 
enables us to derive the term structures of default probabilities (TSDPs hereafter) 
and the recovery rates (RRs hereafter) consistently with rating factors and industry 
factors, where the TSDP is the set of default probabilities that a firm gets defaulted 
by a future time s where s belongs to a positive interval from 0. Our model is a 
cross-sectional empirical model valuing all the given CB prices simultaneously at 
each time. 

In this formulation we take into account the fact that each enterprise has a port- 
folio of multiple business lines corresponding to some different industries, and we 
explicitly incorporate the business portfolio structures of individual firms together 
with rating factors into the model in order to derive the TSDPs and RRs implied 
in CB prices. In our case the portfolio ratios of business lines are measured by the 
sales ratios of industry-wise business lines. Hence in estimation we need the sales 
and rating data together with the CB price data and bond attributes (coupon and 
maturity period). 

Our modeling is based on three fundamental assumptions or viewpoints. The first 
one is about the information content contained in CB prices. The CB prices are 
assumed to be efficiently formed in the market at each time and so reflect or contain 
the investors' views on the term structure of future default probabilities for each 
CB over its maturity period. In fact, their investment decision makings are usually 
based on sufficient information and analysis on firms and hence almost all the CB 
prices should be supposedly consistent with the investors' views. Consequently the 
TSDPs we aim to derive for each firm from all the current CB prices are regarded 
as the investor's forward-looking TSDPs. 

Second, it is viewed in our that the credit condition or credit quality of each 
firm in general depends on its future cash (profit) flows, which in turn depend on 
the business portfolio structure. Therefore, the credit condition or credit quality of 
each firm in general depends on economic trends or business cycles that are often 
different industry- wise. 

The third point is related to the fact that the specification of default correlation 
is very crucial in credit risk analysis. In our CB price modeling, the default correla- 
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tions are naturally introduced from the model structure through considering those 
of stochastic discount functions and the cash flow structure of CBs with defaults. 
On the other hand, in time-continuous setting, the correlations are often assumed 
to be constant though they in fact change constantly and increase in the phase 
of downturn economy where industry-wise business cycles are relevant. Concerning 
modeling the stochastic correlations, in discrete time series analysis Engle (2009) 
extensively treats multivariate return processes with conditional stochastic volatili- 
ties and conditional stochastic correlations. As terminology, when a credit analysis 
contains the credit risks of multiple firms, it is often referred to as multi-name case 
where the correlation structure needs to be specified. Otherwise it is referred to as 
single-name case. 

In applications, the implied TSDP and the implied RR that we derive from the 
formulation of our CB pricing model can be applied to price such credit derivatives 
as CDS (credit default swap). In fact, in this paper we also give a formula for 
valuing a CDS of CB in our discrete-time approach. Since our TSDP depends on 
the business portfolio structure of the issuer, so does the pricing formula of a CDS, 
though the form of the formula itself is rather well known. 

There is a vast literature in the area of credit risk. Most recent researches take a 
time-continuous setting. The books by Duffic and Singleton (2003), Lando (2004), 
and McNeil, Frey and Embrechts (2005) are well known. Most of the articles treat 
a single-name (univariate) case though some recent works consider a multi-name 
case, e.g, Filipovic, Overbeck and Schmidt (2011). But most of the papers do not 
take into account the feedback structure of economies and simply assume exogenous 
processes for credit movements without considering business cycles. Duffie and Kan 
(1996) viewed interest rate processes as dependent and used a state space model 
though they assumed a single Markovian process for the univariate state variable. 
Some papers consider a specific factor for each firm, but the specific factors are 
simply treated as some exogenous processes independently of a common factor. 

In our modeling the stochastic behaviors of the CB prices correspond to those 
of the attribute-dependent discount functions, which are formally expressed with 
attribute- specific forward rates. In association with this viewpoint, Collin-Dufrcsnc 
and Solnik (2001) considers the term structure of default premia in the swap and 
LIBOR markets. Fcldhutter and Lando (2007) decomposes swap rate into a common 
swap rate and swap spreads that are specific to credit or counter party risk, where 
the spreads are referred to as convenience yields. 

The organization of this paper is as follows. In Section 2, some important prob- 
lems are discussed in the formulation of credit risk model to be used in practice. 
In particular, some differences between time-continuous models and our model are 
discussed from a practical view point. In a line with the arguments in Section 2, 
Section 3 presents an extension of the government bond (GB) model that Kariya 
and Tsuda (1994) and Kariya and Tsuda (1996) proposed, where the model is 
shown to have a one-month ahead predictive power for Japanese Government bond 
prices. And we discuss about the relation between the stochastic discount functions 
and spot and forward interest rates which are both attribute-dependent. In Section 
4 the non-defaultable bond model is specified in detail for empirical work and an 
estimation procedure is proposed. The important parts are the specifications of the 
mean discount function and covariance structure of the bond prices, which sepa- 
rates our model from the other models. In this paper GB and non-defaultable bond 
are assumed to be synonymous. 

In Section 5, we formulate the CB pricing model that enables us to derive the 
implied TSDPs and RRs. This formulation is quite different from the models in 
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the literature in that in our modeling the concept of the investors' forward-looking 
views on the TSDPs is introduced and implemented into the model. In Section 6, we 
apply our results to credit risk management in banking and pricing credit default 
swap (CDS) among others. 

2. Problems in credit risk analysis 

In Section 1, concerning credit risk models developed in mathematical finance, we 
discussed about some features of the models from a practical viewpoint. Those fea- 
tures come from the continuous-time setting; Markovian, univariate, instantaneous 
recovery rate, etc. In addition, the following points should be well considered in the 
formulation: 

(1) Data source to be used for credit-risk modeling-past data or current data; 

(2) Default correlation problem-univariate (single-name) model or multivariate 
(multi-name) model; 

(3) Conditional model or unconditional model; and 

(4) Information on industry factors and rating factors. 

The point (1) is important for practical effectiveness of a model so long as the model 
is used for a forward-looking decision-making or investment. In fact, economic, fi- 
nancial and technological environments surrounding firms are evolving constantly 
and rather quickly and hence models using past data on defaults and non-defaults 
and past financial data of firms do not necessarily give us a forward-looking infor- 
mation to get a future credit perspective for each firm. The model based on past 
data may tend to deliver a rather backward-looking information, which may be the 
case of Markovian transition model or hazard rate model, where data on defaults 
and non-defaults over some past period is often used. On this point, it is noted 
that each default event in the past is very firm-specific in its nature, and defaulted 
firms do not exist any longer. Credit risk models that use current market price data 
on existing firms are more forward-looking as the investors try to be rational and 
analytical to make gains. Hence the CB prices that are formed in the market are 
supposed to reflect the views of investors. In this sense our TSDPs whose informa- 
tion is based on current cross-sectional data of CB prices can be regarded as the 
TSDP of investors' perspective views and so it will be more practically effective for 
decision makings. 

Concerning (2) and (3), note that credit risk factors in credit-related instruments 
and derivative products are generally significantly correlated through industry fac- 
tors and business cycles and so credit correlations should be well modeled consis- 
tently in valuation. In addition, business cycles in each industry are often different, 
differently affecting firms that have different portfolios of business lines with some 
overlapping for each other. Furthermore default processes in practice affected by 
such business cycles is in general non-Markovian and mutually dependent. In the 
literature a Merton-type model, which is basically univariate Markovian (condi- 
tional) geometric Brown model for stock (firm value) prices, is often extended to 
a multi-name case where the correlations may be treated as constants with copula 
functions. It is noted that it is not easy to derive the TSDP from stock prices as 
stock does not have a finite time horizon. 

Considering (4), note that the future cash flow (profit) structure of an enterprise 
depends on the portfolio structure of business lines associated with industry fac- 
tors. The cash flow structure matters mostly for the credit quality because most 
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of defaults occur due to the lack of liquidity or cash, not through the imbalances 
in the financial balance sheets of firms. Since the industry factors make a correla- 
tion structure among different credit risks of firms to a great extent, the portfolio 
structure of business lines associated with industries should be taken into account 
for valuing CBs, CDSs and some other credit-related products though the concept 
of industry is relative and needs to be defined in advance. 

Furthermore it is remarked that in reality the recovery rate (RR) after default is 
determined after a long procedure of asset evaluation and negotiation process among 
those of interests, which is costly and very time-consuming. In our formulation it 
is derived as the RR that investors expect at current time and recovery is assumed 
to takes place at the next coupon-paying time, as will be discussed in Section 5. 

3. Pricing non-defaultable bonds 

Suppose that there are G government bonds (GBs) or equivalently non-defaultable 
bonds whose prices are denoted by P g (g = 1, . . . , G). Let t — denote the present 
time and let 

(1) s g i < s g2 < ■ ■ ■ < s gM{g) g = l,...,G 

denote the future time points at which the gth bond generates the cash flows 
(coupons or principal) in view of t = 0. These s g j values are measured in years, 
where s g M( g ) is the maturity period. In this section and the next section future 
time points are measured continuously in years. Assume that the face value is 100 
(yen or dollar) and let c g be its coupon rate (yen or dollar). If coupons are paid 
biannually, the cash flow function C g (s) of the gth bond is expressed as 

{0.5c g s = s gm ,m^ M(g) 

100 + 0.5c 9 s = s gM{g) 
s f s Sgm 

However, in our argument the cash flow function can be arbitrary, so long as the 
future cash flows and their time points are given in advance. 

Let D g (s) be the attribute-dependent stochastic discount function of the gth 
bond defined on < s < s aM ^ with s a M( a ) = max g s gM(p); where the whole 
values of D g (s) are realized all at t = and D g (s) discounts cash flow C g (s) by 
C g (s)D g (s). Under these notations, our basic formulation for modeling G bond 
prices simultaneously at t — is based on the following expression: 

M(g) 

(3) P g = C g (s gj )D g (s gj ) g 1 C. 

3 = 1 

This is an unconditional cross-sectional expression. In (3), we regard the realization 
of price P g as equivalent to the realization of the whole function {D g (s) : < s < 
s gM(g)} with g = 1, . . . , G. Hence the realizations of G bond prices correspond to 
those of D g s am , (m = 1, . . . , s a M( a )] g = 1, • ■ • j G) and the correlation structure of 
these stochastic discount functions implies those of prices. 

In the spot rate approach in mathematical finance a non-defaultable bond price 
at t — is specified as the conditional expectation with attribute-free discount 
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function given the past and present information, which is expressed as follows; for 
each individual price 

M(g) 

(4) P g (l) = £ C a {a gi )b(a gi ) 



with 



D(s gj ) = E 



exp ( — / r u du 



H(r ,s 



93 ' ' 



where {r u : < u < s a M( a )} is a process of instantaneous spot interest rates 
{r u } that is common to all the bonds and Eq[ } denotes the conditional expectation 
given ro at with respect to a risk neutral measure. But the measure is not uniquely 
identified. Here 8 denotes a set of possible parameters when a specific model such as 
CIR (Cox-Ingersoll-Ross) model or Vasicek model is used for the spot rate process 
{r u }. In mathematical finance the conditional expectation can be regarded as being 
taken under a risk neutral measure and (4) is claimed to hold a.s. for all the bonds in 
its no- arbitrary theory. In fact, conditioning variable ro is the only random variable 
making all the G bond prices realized. Hence in reality it does not follow that (4) 
holds a.s. as the bond prices are not functionally related and in fact maturities 
and coupon rates are different. In such a spot rate approach, modeling the spot 
rate process yields the conditional discount function through which the zero yield 
curve {Ru : < u < s a M( a }} defined by H(ro,s,9) — exp(— R s s) or equivalently 
R s = (— 1/s) log H(r , s,9) is obtained. In the sequel we use the real measure that 
generates real data. 

Another remark on this approach is that this specification of spot rate process for 
bond-pricing ignores such bond attributes as coupon rate or maturity. Empirically 
speaking, it is often observed that bond prices formed in the market depend on 
such attributes, in which case it is required to take into account such attribute- 
dependency in the specification of the spot rate process. In modeling a swap rate 
process Collin-Dufresne and Solnik (2001) and Feldhuttcr and Lando (2007) take the 
dependence of swap rates on credit attributes into account and specify a swap rate 
process as the sum of an abstract risk-free rate process {xi s } and a convenience yield 
process {x2 gs } where they are assumed to be independent. Here the convenience 
yield represents such attributes as liquidity premium, credit premium (collateral 
condition), etc.; 

^gs ^ls ~t~ ^2gs • 

The attribute-dependent convenience process {x2 gs } can play an adjusting factor 
for fitting the model as it can be arbitrarily specified. Using this process the discount 
function in (4) becomes attribute-dependent; 

exp - / r gs ds 



D g {s gj ) = E 

In this paper we do not use this spot rate approach in (4) with conditional expec- 
tation but take a forward rate approach with unconditional expression, and make 
it the following attribute-dependent model as in (3); 



(5) P g = 2^ Cg{sgj)D g {sgj) with D g (s) = exp ^- f gu du^ 
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where {f gs : < s < s a M( a )} is an instantaneous forward rate term structure 
process whose values are realized all at t = 0. In other words, for each g a realization 
of P g corresponds to that of the whole path {f gs : < s < s a M(a)}- Note that (5) 
is equivalent to (3). 

On the other hand, as an attribute-free forward rate process one may use the 
time-continuous HJM (Heath, Jarrow and Morton (1992)) model, which describes 
for each individual g a process of term structures {ft s '■ < s < St a j\/(a) , t > 0} 
with the discount function Dl(s g j) — exp (— J^ 93 ft s ds) attribute-free. The HJM 
model is specified conditionally and the Markovian expression is usually explored 
with the no-arbitrage argument. Even in this case, one single path {f s : < s < 
SaM(a)} does not make (5) hold a.s. for all g either and so we use (5). 

Finally it is remarked that a specific cash flow pattern guaranteed by holding 
some GBs will be a big value to such institutional investors as pension funds or 
life insurance companies because they need to match cash inflow with cash outflow 
over a long time horizon. In other words, coupon and maturity are important at- 
tributes which affect investment decisions with a future perspective. For example, 
depending on cash inflow-outflow structures of investors and on future perspectives 
on movements of interest rates, it may happen that a GB of 2 year maturity and 5% 
coupon is less preferred to a GB of 6 year maturity and 3% coupon. It is noted that 
those institutional investors do not necessarily prefer de-coupon or stripped bonds 
engineered by investment banking as making a portfolio from these stripped bonds 
to match the cash inflows and outflows is costly and involves additional credit risk 
of investment bankers. 



4. GB pricing model 

In this section, we implement an attribute-dependent GB pricing model with a 
stochastic discount function D g (s) in (3) or (5). First, let D g {s) be decomposed 
into the mean function and the stochastic deviation function as 

(6) D g (s) = D g (s) + A g (s) . 
Substituting this into (5), it follows that 

M(g) M(g) 

(7) Pg = ^ Cg(S gm )Dg(S gm ) +T]g T) g = C g Ag = ^ Cg ( S gm ) Ag ( Sg m ) , 

m—1 m— 1 

where 

(8) 

Cg = (Cg{S gl ), Cg{Sg M (g)))' '<^d Ag = (Ag(S gl ), Ag(S gM{g) ))' . 

The expression (7) corresponds to the case in (5), but without specifying the 
attribute-dependent forward rate process {f gs : < s < s a M(a)}j the- correspond- 
ing mean discount function D g (s) on [0,s a M( a )] is assumed to be continuous in s 
and then it is uniformly approximated by a pth order polynomial; 

(9) D g (s) = 1 + (Sixzig + 5 12 z 2 g + S u z 3g )s H h (S pl zi g + S p2 z 2 g + S p3 z 3 g)s p , 

where Z\ = 1, z% = c g , and z 3 = s g M( g ) are the attribute variables of the gth bond. 
In this specification the parameters are common to all the mean discount functions 
for g = 1, . . . , G and hence they are estimable with G bond prices, so long as G 
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is greater than the number of the parameters contained. Substituting (9) into (7) 
yields 

M(g) 

Cg{Sg m )Dg{s gm ) = a g +5\ 1 dg\ \ +<5;L2^gl2 +<5l3^gl3 +" ' ■+5 p ldg p i+S p 2dg p 2+5 p 3d gp 3 . 

m— 1 

where 

M(g) M(g) 

a g = ^ t C"g( s gm) and d g ij = ^ ] Cg{ s gm) z gjS gm ■ 

m— 1 m— 1 

Here « in d 9 ,-j denotes the attribute suffix and j the polynomial order. Thus letting 

Xg = {dgll, dg21, dg^\] dg\2, d g 22, dg32', • • • ; dg\p, d g2 p , dg^p)' 

and 

X = (xi,x 2 , . . . ,xa)' 

we have a regression model 

(10) y = Xp + Vl 

where y = (y x , y 2 , . ■ ■ , ye)' with y g = P g - a g , r} = (%, . • . , VgY and 
/3 = (<5ii, 5i 2l <5i3 : (52i, ^22, ^23; ■ • • > 5 P i,3 P 2, ^3) . 

In (10) the specification of the covariance matrix rj is crucial since specifying the 
covariance structure of P = P\, . . . , Pg)' or equivalently the covariance structure 
of -q stochastically describes a structure of the joint realizations of G bond prices. 
In view of (7) the specification is directly related to that of the covariances of the 
stochastic discount factors D g (s g j) and Dh(sh m ) at each cash flow point s g j and 
Shm of the gth and hth bonds. We specify it as 

(11) COV(Dg(Sgj), D h (s hm )) = <J 2 Xg h fg h . jm , 

where a 2 is a common covariance and covariance factor, \ g h is a covariance part 
related to the differences of maturities and f g h-jm is another covariance part related 
to the difference of the cash flow points s g j and Shm- These two parts are further 
specified as 



(12) Xg h = 

with 



pe gh g ^ h 
e gh = exp (-£\s g M(g) - ShM(h)\) 



and 

(13) fgh-jm = exp (-6\Sgj - s hm \) , 

where we assume that < 9, p, £ < 1. These specifications imply; 

(1) as is expressed in e gg of X gg , the longer the maturity of each bond is, the 
larger the variance of each price is, 
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(2) as is expressed in e g h of X g h , the larger the difference of the maturities of two 
bonds, the smaller the covariance is, and 

(3) as is expressed in f g h-jm, the closer the two cash flow points are, the larger 
the covariance of the discount factors D g {s g j) and D^Shm) is. 

Under this specification, the covariance matrix of rj is given by 

(14) Cov{ v ) = (Cov(r) g ,r) h )) = (Cov(P g ,P h )) = a 2 (X gh cp gh ) = <r 2 $(0,p,O 
with 

M(g) M(m) 

Vgh = Cg{ s gj)Ch( s hm) fgh-jm ■ 

j—1 m= 1 

As in Kariya and Kurata (2004), the unknown parameters are efficiently esti- 
mated by the GLS (generalized least squares) method, in which we minimize 

(15) e, P ) = [y- xp}' me, P , or 1 [y - xp\ 

with respect to the unknown parameters. First, for given (#,/},£), the minimizer of 
this function with respect to p is known to be the GLSE; 

p(6, P ,0 = [x^ie^p^r'xr'x^ie^p^y'y 

and then the marginally minimized function ip{P, 9, p, £) with substitution P(Q, p, £) 
is minimized with respect to (0,p,£), yielding the GLSE (f3,9,p, £) where a grid 
point method for split points of (9, p, £) may be used. 

In Kariya and Kurata (2004, pp 55-63), an empirical performance due to Kariya 
and Tsuda (1994) is demonstrated as an example of GLS estimation where the order 
of the polynomial in (9) is set to 2 with z\ = 0, z-i = coupon, and zj, = maturity. 
Even in this case, the residual standard deviations are 0.338 yen for December 
27, 1989 with G = 70 and 0.312 yen for January 31, 1990 with G = 70 where 
the face value is 100. Recently, the effectiveness of this model is empirically and 
comprehensively tested with Japanese bond data in Kariya et al. (2011). 

It is noted that in the specification of the attribute-dependent mean discount 
function in (9), setting z\ = 1, z% = 0, and Z3 = yields the attribute-free specifi- 
cation. In this case the parameters involved are estimated in the same manner to 
get the attribute-independent mean discount function D(s). This D(s) is converted 
to a yield curve by R s = — s _1 log D(s), which is often referred to as a risk-free 
yield curve. This attribute-free discount function may be used to price a CDS in 
Section 6 because the cash flows in CDS does not reflect the same cash flow pattern 
as the corresponding CB. 

5. CB pricing model for deriving TSDBs 

In this section we propose a formulation of the CB pricing model that enables us to 
derive the TSDPs and RRs for each pair of industry index and rating index where 
(1) CB price data, (2) data on industry-wise sales ratios of firms that issued the 
CBs and (3) credit rating data are assumed to be given at t = 0. As has been stated 
in Sections 1 and 2, in our modeling the differences of the business line portfolios 
of firms are taken into account in terms of sales ratios of the firms. As a matter of 
a fact, each firm has different exposures to industry- wise business factors and so it 
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has a different profit structure in association with industry business cycles, which 
is greatly relevant to analyzing credit qualities of CBs. 

To formulate our model, suppose that at t = there are K CBs to analyze 
and let {s M ; I = 1, . . . , M(k)}, k = 1,2,..., K with s kl < s k2 < ■ ■ ■ < s k u(k) 
denote the future cash flow time points of those CBs as in the case of GBs. Also 
let s Q M( a ) — max k SkM(k) an d let the cash flow function C k (s) of the fcth CB be 
defined on < s < s a M( a ) for all k though it is zero except for the above finite 
points. 

On the other hand, if the firm that issued a CB gets defaulted before its maturity 
the coupons to be paid thereafter are not paid and some portion of the face value 
100 may be paid after a long procedure of legal and practical settlements, where 
the portion relative to 100 is called recovery rate (RR), which is in general of a 
stochastic nature. The expected or empirically averaged value of the RR is known 
to depend on its credit grade via a rating agency. Hence actual cash flows from 
a CB depend on how likely a firm that issued a CB is to get defaulted and what 
the expected RR is when it gets defaulted. In the market the TSDPs are evaluated 
simultaneously and consistently together with the expected RRs and are implicitly 
reflected in their market prices of CBs. 

5.1. Basic formulation of CB pricing model 

On this viewpoint, let T k be the first random time of the fcth CB to default and let 



Then L ks defined at t — is the indicator function of default event {r k < s} and 
identifies if the kth CB gets defaulted before or on a future time point s. Then 
the actual cash flow function C k (s k j) at a future coupon generating time s k j is 
expressed as 



This means that C k (s k j) — C k (s k j) if L kSkj = 0, or equivalently if the firm has 
not defaulted until s kj , and C k (s kj ) = 100^(i(k)) if L kSkj {\ - L kSk] _ t ) = 1, or 
equivalently if it had not defaulted at s k j-i and has defaulted at s k j. Here 7(z(fc)) 
is the mean RR of the kth bond with rating grade i(k), where the credit rating is 
indexed by natural number i = 1, 2, . . . , / and the smaller the number is, the higher 
the credit grade is. Here this specification implicitly assumes that the recovery 
payment is made at s k j if a default event occurs in interval (s k j-±, s k j]. Note that 
s k j — s k j_i is typically about a half year. 

However, the expression (17) itself exhibits the future relation and has never been 
realized at t = 0. Therefore investors expected cash flow at t = is formulated as 
for interval (sjy_i, sjy] 



(16) 




(17) 



C k (s kj ) = C k {s kj )(l - L kSh .) + 100 7 (z(/c))L fcSfcj (l - L kSkj _J . 



(18) 



Ck(skij) = C k {s k ,j)[l -p k {s jk : i(k))) 

+ I00~f(i(k))[p k (s kj : i(k)) -p k (s kj ^i : i(k))]xk(skj) , 



where p k (sj k : i(k)) — E[L kSkj ] is the default probability corresponding to the 
event {r k < s k j}, which is the probability that the fcth CB gets defaulted before or 
on s k j and Xk(skj) is the indicator function of the cash flow time points {s k f,l = 
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1, . . . , M(k)} of the fcth CB. With investors expected cash flows in (18) we formulate 
our CB model as 



(19) 



M(k) 

Vfc = ^ Ck(skj)D k (skj) 



Before we proceed further, it is remarked that a typical model in the mathemat- 
ical finance literature is the conditional univariate model that introduces the credit 
element into the discount part; 



Vu 



M(k) 

C k (s k3 )E 



3=1 



exp 



(r u + \ k u )du 



where {r u } is a Markovian spot rate process and {Ajj} is an instantaneous Marko- 
vian default intensity process that discounts cash flows together with the spot rate 
process. In this expression it is assumed that {r u } is independent of {L ks } , which 
is unlikely but yields 



En 



(1 



)exp M 



. du 



= En 



exp 



(r u + A*) du 



by the Doob-Meyer Theorem and hence the default intensity {A„} are assumed to 
be exogenously independent of {r u } . It is often the case that such processes as CIR 
model are assumed for {A u } , the discount function E [] is evaluated analytically 
and then the unknown parameters therein are calibrated or estimated with the 
present or past data of the fcth CB only. This model is sometimes extended to a 
multivariate case where the correlation of A^ and A^ is often assumed to be constant. 
Clearly this modeling approach is quite different from ours. 

Another remark is that the expression (18) may be regarded as the conditional 
expectations of investors and in that way the discount function in (19) and hence 
the model (19) itself may be regarded as a conditional expression to get a dynamic 
model. Now coming back to our case, the expression in (19) corresponds to the one 
in (5) and hence the rest of the argument is similar to the non-defaultable case. 
That is, the stochastic discount function is decomposed as 



(20) 



D k (s) = D k (s) + A k (s) 



and as D k (s) we use the mean discount function estimated with G government bond 
prices where the attributes of the /cth bond are inserted into the discount function 
in evaluation. Consequently 



(21) 



{D k (s kj ) : j = 1, . . . , s kM{k) ■ k=l,...,K} 



is a set of known values and it follows from (19) and (20) that the CB pricing model 
is 



(22) 



V k = 



M(k) 

E 

3 = 1 



C k (s kj )D k (s k j) +E k with e k 



M(k) 

E 

3=1 



C k (s kj )A k (s kj ) 



It is remarked that a joint model of GBs and CBs can be formulated in which the 
common mean discount functions are estimated simultaneously, though we do not 
follow this because of its complexity. Note that there are about 3,000 CBs in the 
Japanese market and about 50,000 CBs in the US market. 
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5.2. Specification of TSDPs with business portfolio structures and 
credit discounts 

Next we specify the default probability function of the kth bond in (18) as 

J 

(23) P k (s : i(k)) = Y, MjMs ■ i(k),j) , 
where w k {j) > 0, w k(j) = 1. Here 

(24) P(s:i,j) i = l,...,I, j = l,...,J 

is the generic or common TSDP with credit grade % and industry j , which is 
independent of specific CBs. In this paper it is assumed to be approximated by a 
polynomial of the qth order; 

(25) p(s : i, j) = a\ J s + a^s 2 + • • ■ + ap q . 

On the other hand, {tUfc(l), . . . , Wk{ J)} in (23) is the set of the sales ratios of the kth 
CB issuer corresponding to the industry indices j = 1, . . . , J and it is regarded as a 
business portfolio in terms of industry-wise sales, where the industry categorization 
is determined in advance. 

Now from (23) and (25) the TSDP of the fcth bond is expressed as 

(26) Pk {s : i(k)) = sw' k a\ + s 2 w' k a l 2 +■■■ + s««/ fc a*' 

= {sw' k ,s 2 w' k , . . . ,s q w' k )f3(i) 
= w k ( S )'(3(i) 

where 

Wk = (ttffe(l), • ■ . , Wfe(J))', 
%• ( il i.J\f t -i 

a h = \ a h >■ ■ ■ , a h ) h = l,...q, 

w k (sY = (sw' k , . . .,s q w k ) . 
Also the expected cash flow function in (18) is expressed as 

(27) C k {s kj ) = Ck(skj) + Zk(skj,skj-i ■ y(i(k)))'P(i) , 
where 

z k (skj,s k j-i : l(i(k)))' = -C fe (s fej -)w(s fci )'+1007(i(fc))xfc(sfc : ,-)[ ti; ( s fei)'^ w ( s fcj-i)'] 
Therefore the pricing model in (22) is given by 

Af(fe) M(fe) 

V k = Y Ck(s k j)Dk(skj) + Y D k (s kj )z' k (s kj ,s kj -i : y(i(k)))P((i(k)) + e k 

3=1 J'=l 

= P k + [uk + y(i(k))v k ]' + e k , 
which yields a regression model 

(28) yf ] = [u k + y(i{k))v k ]'j3(i) + e k 



Corporate bond model 



13 



with 

V k -P k , 

M(k) 

Ck(skj)D k (s k j) , 

3=1 
Af(fc) 

- C k{ s k3)Dk{s kj )w k {s kj ) and 

3=1 
Af(fc) 

100 D k (s kj )xk(skj)[wk(skj) ~ Wfe(s fe j_i)] . 

3=1 

In this expression, is regarded as an expected (or a theoretical) GB price with 
non- default able cash flow {C k (s k j) ■ j = 1, • • • , M(k)}, and hence y£® is the 
difference between the fcth CB price with credit grade i{k) and the corresponding 
non-defaultable bond price. The difference yl tends to be non-positive, because 
C k (s k j) < C k (s k j) and a CB is of the less creditability and less liquidity than the 
corresponding GB. We simply call y fe the credit risk discount of the fcth CB price. 

Consequently (28) forms a regression model for the credit risk discounts. In fact, 
supposing that there are 7Q CBs of credit grade i and letting y(i) = (y\, . . . , y l K . )' , 

Xi(i) = (u\ , . . . , u % K .)' and X 2 (i) — (v\ , . . . , v % K . )' the pricing model for CBs of 
credit grade i is reduced to a regression model; 

(29) y(i)=X(i, 7 (i)))8(i)+e(i) 

with X(i,~f(i)) = X\{i) + j(i)X2(i), i — 1, .. . ,7. Here K{ > 2Jq is assumed for 
the identifiability and estimability of /3(i) and 7(1). In fact, Ki > 2Jq is a neces- 
sary condition for the uniqueness of (/?(?), 7(1)) in the sense that X{i,^j{i))P{i) — 
X(i,"/{i)*)P(i)* implies j3{i) = and = It follows from (29) that the 

credit risk discount vector y(i) of the ith credit grade is explained by the regres- 
sion matrix X("/(i)) which depends on the unknown RR ~f(i). Combining all the 
regression models over i = 1, . . . , 7 yields 

(30) y = Xfi + e 

with X = [X(i,j(i))] and [X(i,j(i))] is the K x IJq block-diagonal matrix with 
the ith block matrix X(i, and /3 = (/3(1)', . . . , /3(7)')' with K = K 1 + --- + K I . 

5.3. Specification of covariance structures 

To specify the covariance structure, we write the dependency of the mean cash flow 
function on unknown parameters as C k (s k j) = C'k{s k j ■ 7(*(&)))- Then the 
error term in (22) is given by 

M(fc) 

e k = °k(s kj ■ P(i),7(i{k)))A k (s kj ) . 
3=1 

Hence the covariance of two error terms of the fcth and ith CB prices with rating 
i(k) and i(l) respectively is assumed to be 



i(k) 
Vk = 

Pk = 

Uk = 
Vk = 
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(31) 

M(k) M(m) 

Cov{e k ,ei) = X! Ck(s k j ■ /?(*)> l(i(k)))Ci(s lm : f3(i), j(i(l)))Cov(A kj , A,. 

.7 = 1 m—l 



where 

with 
and 



^fei = ( Pueki k^l, i(k) = = i 

^Pij&ki k^l, i(k) = i, = j, i ^ j 

e M = cxp (-^ |s feM (fe) ~ s iM(i)\) i(k) = i, = j, 



M(k) M(m) 
j — 1 m—l 



with 



&fci-jm = cxp (-0 |s/y - S; m |) . 

In this specification, when the fcth CB and the Ith CB are of the same credit rate, 
i.e., i{k) — i(l) = i , then the covariance is of the same form as the one in the non- 
defaultable case. If the two CBs are not in the same credit category, ptj will account 
for the cross-correlation between the two categories. Using the above specifications, 
the covariance matrices for regression models are obtained; 

Cov(e(i)) = (Cov(e k ,ei)) = <T 2 (\ k Wki) = cr 2 $(/3(i), pa, £u) = cr 2 $ ri 

(32) 

Cov(s(i),e(j)) = CT 2 *()8(*),/80').7(*).70').P«.f«) = * 3 and 

Cov(e) = (Cov(e(i),s(j))) = a 2 {^ 3 } . 

(33) 

In the CB case, the covariance matrices depend on the regression coefficient f3(i)'s. 
5-4- Estimation procedure 

Under these formulations we propose the following grid procedure to simplify the 
model estimation. 

(1) For each credit category, the GLS estimation is pursued: fix i. 

(a) For each given (7(4), pa, £u), each of which moves over 0, 0.1, . . . , 0.9, 
f3(i) is estimated by a repeated procedure. Setting (5(i) = in <&a and 
minimizing 

i, = [y(i) - X(i, 7 (i))/3(i)]'[<M/3«, 7(<). Pa, tuT'lv (<) - 7«)/3«] 
yields the first step GLSE with <I>-° ) = $u(0, pa, 

= ?( ! ,7i'*f 1 i(i,i( ! ))r 1 i(M( ! ))'$f 1 !/(i) • 
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Substituting this into to get $|p = ^(^(i)^, 7 (i), pu^a), 

applying the same procedure yields the second step MLE f3(i)( 2 \ Re- 
peating this procedure a couple of times with a performance evaluation 
rule, we obtain an approximate minimizer given 

(b) Repeating the procedure in (a) over possible values of (7(2), pa,!;™)-, we 
obtain (/?(?)*, 5?^, 7 (i)*, p*j, that minimizes the objective function ap- 
proximately. Once we obtain this, we may repeat (a) in the neighborhood 
of the approximate optimal value by splitting the neighborhood with a 
finer division. 

(2) A simultaneous estimation procedure for the combined model is obtained as 
follows. Fix (P(i)*,j(i)*,p*i,i*i) i = 1, • • • ,1 in the matrix 

®ij = *G8(*).^0')i7(*).7(j)./»tfiftf) Ji 

meaning that the full covariance matrix in (32) carries the /(/ — 1) unknown 
parameters {pij,£ij ■ i j}- Hence in the same way as in (1) we minimize 

to obtain the first approximate simultaneous GLSE (3^u- One ma y replace 

/3(i)*'s in the fixed 0(i)*, p* u , |* ) % = 1, . . . ,7 by /3^ H and repeat the 
minimization procedures again to get the second approximate simultaneous 
GLSE fifliv 

(3) This minimization procedures are repeated over the degrees q = 2, ... ,6 of 
the polynomials in (25) to get an optimal simultaneous GLSE Pf u u- 

With this model, an extensive empirical research is under a way with Japanese CB 
data. A simple empirical work is demonstrated below in a limited model. In this 
model we do not take into account the fact that each firm has a specific business 
portfolio over industries. In other words, we set J = 1 and pick the AA class of 
208 CBs with the assumption of recovery rate 7 = 0. A plot (see Figure 1) of the 
residuals {i^ A = Vk — Vk} shows that even in this simplified model the performance 
are rather good and hence our full model will be more promising as a model to derive 
the TSDPs. In fact, the standard deviation of et A 's is 0.419 yen. 
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Plot of residuals for AA: no industry distinction, y=0 
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6. Applications 

Once a full estimator f3f u u is obtained, the following applications are considered. 

[1] The whole implied TSDPs in (25) for each pair of credit grade index and 
industry index are given by 

(34) p(s : i,j) = c^s + a^s 2 + ■■■+ 6$ s 9 i=l,...,I j = 1, J 

together with the implied RRs ■ i = 1, • • . ,1}- These TSDPs and RRs are 

basic information sources for credit-related products. 

[2] The TSDP of the fcth CB or its issuer with credit grade i(k) and business 
portfolio {w k (l), . . . ,Wk(J)} over J industries is obtained from [1] as 

,/ 

(35) p k (s : i(k)) = ^2w k (j)p(s : i(k),j) . 

3 = 1 

This is the TSDP {P{r k < s) : s > 0} Of the issuer of the fcth CB where r k is the 
default time as in (16). The expected RR is 7(i(/c))* = 7(z) with i(k) = i. 

[3] Inserting (3f u u into (28) and using (27), the credit risk discount level of the 
fcth CB or its issuer with credit grade «(fc) and business portfolio {u>fc(l), . . . , Wk(J)} 
over J industries is given by 

[v>k + j(i(k))*v k ]'j3fuu 

M(k) 

^ D k (skj)W(p k (skj ■ i(k)),Pk(skj-i ■ i(k));j(i(k))) 

3=1 

with 

(37) W(s kj ) = W{ Pk {s k3 : i{k)),p h {s hj . x : i(k));j(i(k))) 

= [100j(i(k))xk(skj) ~ C k (s kj )]p k (s kj : i(k)) - 100 7 (i(fc)) 
x Pk(sk 3 -i ■ i(k))xk{skj) 

where the discount is relative to the fair value P k = V k + (— j/L ) of the correspond- 
ing non- default able bond. Here P k is the model value of non-defaultable bond with 
the same coupon and same maturity as those of the fcth CB. Consequently the fair 
spread rate of the fcth CB relative to the corresponding (attribute-dependent) GB 
is given by 

(38) s Pk = (P k - V k )/P k = -y k (k) /P k ■ 

Hence if the actual spread (P k — V k )/P k = —y k ^/P k is too low or too high relative 
to this fair spread, one may use the information for investment decision making 
with the expectation that the actual spread eventually converges to the fair level. 
[4] CB Portfolio Management. The fair value of a portfolio {tu k } of CBs is given 

by 

K 

(39) ^Vk 
fe=i 



(36) yf > = 



K M (a) K 

— J^a7fc C k (s am )D k (s am ) = ^n7fc[P fe + y k ] 

k =l m=l fe=l 
M(a) M(a) 
= ^ [^-am + B arn ] — S C a 



-'am 

m—1 m=l 
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fc=l 



K 



where {tUk} is a portfolio of CBs in terms of physical units, 

K 

(40) A am = 
and 

(41) B am = ant 

)W(p k (s am : i(k)) 

i Pk \ s am—l 

fc=i 

Here A am is the present value of the cash flows at future time s am if none of the 
bonds are defaulted before or on s am where {s am } is the combined set of cash flow 
points, while B arn is the present value of the expected loss at s am , which is generally 
negative. 

What is important in this expression is that these are regarded as fair present 
values of those of the corresponding loan portfolio in a bank since a fixed-rate loan 
is regarded as a CB and a variable-rate loan may be converted to an equivalent fixed 
loan via interest swap rate at the time of evaluation. Though we do not elaborate it, 
based on this expression, the portfolio can be adjusted to a desirable loss structure 
relative to industries and rating grades via 

J 

Pk(s;i) = 2^Wfc(j)p(s : , 

3=1 

where (37), (39) and (41) are combined to sort out the overlapping structure of 
firms portfolios and extract the industry positions for each p(s : 

Another example is about the durations of expected losses. Putting b m = B am / B 
where 

M(a) 

B = 2_j Bam , 

m—1 

{bm} represents the relative sizes of losses distributed over {s am } and hence 



is the average duration of the losses. On the other hand, {a m } with a m — A am /A 
and 

M 

a = ' A am , 

represents the relative sizes of the cash inflows distributed over {s am } if none of 
the CBs or loans get defaulted, and so the duration is 



dm Sam 



E 



The two durations are in general different and the expected actual duration becomes 



E 



m J am 
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with Cm = C am jC and 

M 

C = ^ ' Cam j 
m— 1 

The above argument can be extended to the case where the portfolio contains GBs 
as well as CBs. 

[5] CDS Pricing. Now we apply the result to the valuation of credit default swap 
(CDS). The CDS is notorious as a credit derivative for enlarging and broadening 
the financial crisis because the protection sellers like AIG sold CDSs more than 
their capacities. However CDS itself is an important instrument because it is an 
insurance-type product (contract) guaranteeing the principal of a specific CB (or 
loan) for the protection buyers by covering the loss when the issuer C of the CB gets 
defaulted. A protection (insurance) buyer B in this contract is the one who pays the 
protection premiums quarterly or biannually to a protection seller A over a certain 
period till its contract period or default of C in the period in order to protect the 
principal. The CDSs of this type are quoted in the market for trades of credit risk 
on GBs and CBs. The quotation includes the ask-side as well as the offer-side of 
protection against each bond, from which a market value of the credit risk of a 
CB is found. Here we first give a price (premium) formula for CDS in a discrete 
time setting, which is different from those in a continuous setting. As is stated in 
Section f, a discrete time setting is more suitable in credit risk analysis than a 
continuous time setting. For example, a credit variation process is not Markovian 
and a settlement process after a default takes often 3 months through 6 months. In 
the discrete-time setting, the uncertainty of the time delay for settlement is directly 
incorporated into the model if necessary. 

Now to price a CDS, let m = 0, 1, ... M be daily counted time points at where 
the time unit h of day is measured in year as h = 1/365. Time point to corresponds 
to mh years from 0. Let 

(42) T = {0 < mi < to 2 < • • • < m k } 

be the biannual premium payment time points counted in days viewed at 0, and 
at these time points the protection buyer pays premium x for principal value 100. 
Further let N c be the default time of firm C as in (16) and let the default process 
be denoted by 

{L m } with L m = X{NC< m } ■ 

To derive a theoretical premium x , we distinguish the two cases; to T or to € T. 

When to T the cash flow at m of the protection buyer B occurs from the 
protection seller A only if firm C gets defaulted at to. Hence the payoff of A is 
expressed as 

(43) Vi = -100(1 - L m -i)L m + 100 7m (l - L TO _i)L m . 

In (43), if C had defaulted at to — 1, L m -\ = 1 and so = , i.e., no cash flow 
at to. Here this expression assumes that on the day of default, protection seller A 
pays 100 and receives recovery 1007 m , which is not realistic, but we can change it 
in a way as we wish. In fact, if the protection seller A pays the principal to B in 
14 days after default and A gets the recovery money in 90 days after default, the 
expression can be modified as for the payoff of A at to 
(44) 

U m = -100(1 - L OT _i_i4)i m _i4 and U m+90 = 1007„ l+9 o(l - L m _i_i 4 )L r „_i 4 . 
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However, for simplicity we here assume (43) in which the payments are made at m 
when the default occurs in (m — l,m]. 

When m 6 T, B pays lOCte to A if the default has not occurred until m and 
hence combining this with (43) the payoff of A becomes 

(45) U A = -100(1 - 7m )(l - L m - X )L m + 100s(l - L m ) . 

Note that by the definition of L m , for m < n, L m — 1 implies L n — 1 and hence 
U A = 0. Now to value a theoretical premium a; in a way as in no-arbitrage theory, 
each payoff U m at m is regarded as a derivative and is valued relative to cash 
managed up to m with daily spot rate process {rj}; 



(46) B m = exp 

Then the martingale condition gives the value at of U n by 

(47) v A {m) = E* [U m /B m ] = E* [d(m)U m ] , 

where in the theory the measure is risk neutral though it is not uniquely identifiable. 
If the interest rate process is independent of U n , which would not hold in reality, 
the value is expressed as 

(48) v A {m) = 100L»(m)[-(l - j(i(C)))Q(N c = m) + 5{m)xQ{N c > m)] , 

where S(m) = 1 if m G T, and it is otherwise. Here we use an attribute- free 
discount function; 

D(m) = E [l/B m ] = E [d(m)} . 

This discount function can be chosen as the one stated at the end of Section 4. 
Therefore the expected payoff at of A is 




m— 1 

On the other hand, the payoff of B should be fair to that of A, meaning that it is 
equal to V A and V B = — V A . Consequently the premium x must satisfy V A = 0. 
The CB-CDS premium described above is valued as a fair value; 

(49) x = (i-7WC)))£:j(HQ(^ c ^) 

Ef=iQ(^ C > m k )D(m k ) 

where 7(1(6")) is the recovery rate of the CB issuer C and {mk} is the set of premium 
payment times. In our case, j(i(C)) is the implied recovery rate and the default 
probability is the implied TSDP, 

,/ 

Q(N C < m)=p c {mh : i{C)) = Y J MC)p{mh : i(C),j) , 

3=1 

where the industry-wise business structure of the issuer C is taken into account. In 
(36), one may replace 

M 



Q{N C = m)D c {m) 
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by the continuous time expression 

Sm(C) 

D(s)p' c (s:i(C))ds, 
where p' c (s : i{C)) is the derivative of the TSDP pc- 
7. Conclusion 

In this paper, we formulated a CB pricing model from which the TSDPs are derived 
for each industry index and each rating index and the recovery rates are derived for 
each rating index. A notable feature of this model is to take into account the fact 
that each firm has a different portfolio structure of business lines over industries. 
The TSDPs and recovery rates are basic inputs in credit derivatives and credit risk 
management. Though we wait for a thorough empirical work, the model itself and 
derived outputs will be very useful for practical decision makings. Further these 
results can be also applied to pricing some multi-name credit derivatives. 
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